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I compute from first principles the local heating rate q (the amount of electromagnetic energy converted to 
heat per unit time per unit volume) for electromagnetic waves propagating in magnetically and electrically 
polarizable media. I find that, in magnetic media, this rate has two separate contributions, q( v > and q^ s \ the 
first coming from the volume of the medium and the second from its surface. I argue that the second law of 
thermodynamics requires that the volume contribution be positive and that this requirement, in turn, prohibits 
negative refraction. This result holds for active or passive media and in the presence of anisotropy and spatial 
dispersion. © 2008 Optical Society of America 



1. Introduction 

Macroscopic electromagnetic theory of material media 
which can simultaneously support electric and magnetic 
polarizations denoted by P and M, respectively, has 
been developed over a century ago and is exposed in 
many standard textbooks. However, in the optical fre- 
quency range and at higher frequencies, this theory has 
long been viewed as purely abstract and noncmpirical. 
Even at much lower frequencies, materials which simul- 
taneously exhibit nonzero magnetic and electric suscep- 
tibilities (and are sufficiently transparent to allow any 
noticeable penetration of electromagnetic field into their 
interior) are quite rare and exotic. 

While it is possible to argue about the physical attain- 
ability of artificial materials with nonzero electric and 
magnetic susceptibilities in any given frequency range, 
nothing precludes us from formally developing the elec- 
trodynamics of such media based on the macroscopic 
Maxwell equations. In particular, this approach was 
adopted by Veselago in the now famous paper (Ref. 1). 
Veselago was interested in materials whose electric per- 
mittivity e and magnetic permeability \x are simultane- 
ously negative and which can exhibit the so-called neg- 
ative refraction - the physical effect which takes place 
when an electromagnetic wave entering the medium 
(e.g., from vacuum) is refracted at the "negative" Snell's 
angle. 

After the publication of a more recent paper by Pendry 
in which a perfect (subwavelength-focusing) lens built 
from a negatively-refracting material was proposed [2], 
enormous attention was attracted to negative refraction. 
Numerous proposals for manufacturing artificial materi- 
als with negative refraction have been put forth. There 
has also been a rigorous effort to demonstrate negative 
refraction experimentally; see, for example, Refs. 3-5 and 
references therein. 

Simultaneously with the activities mentioned above, 
there has also been a persistent effort to subject the 



physical attainability of negative refraction to doubt. 
Perhaps, the most consequential of such exploits is the 
recent paper by Stockman [6] in which it is shown from 
the causality principle that, in a negatively- refracting 
material, the rate of dissipation of the electromagnetic 
energy into heat can not be lower than a certain thresh- 
old and that any attempt to compensate for such dissi- 
pation, e.g., by introducing optical gain, will necessarily 
destroy the negative refraction. It is worthwhile to note 
that low dissipative losses are essential for the realization 
of the original Pcndry's proposal for the perfect lens. 

In this article I confront the phenomenon of negative 
refraction with another fundamental physical principle, 
the second law of thermodynamics. I show that the gen- 
eral requirement for negative refraction in isotropic me- 
dia, namely, 



Im(e^) < 



(1) 



is in contradiction with the latter. To do so, I compute 
the heating rate q(r) in a magnetically and electrically 
polarizable medium. I do this by two methods, one in- 
volving the expression — V • S, where S is the Poynting 
vector, and the other involving the expression J- E, where 
E is the electric field and 



dP 

~dt 



cV x M 



(2) 



is the total current induced in the medium (we assume 
that there are no external currents or charges). Quite un- 
expectedly, I obtain different results. I claim that the ex- 
planation for this discrepancy is that the Poynting vector 
in a magnetically polarizable medium must be defined by 



S = — E x B 

47T 

rather than by the commonly used formula 
_ c 
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E x H 



(3) 



(4) 
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Arguments for the validity of ([3]) are given below. 

When the definition ([3]) is adopted, the two methods 
of computing q(r) give the same result. It further turns 
out that the heating rate has two separate contributions: 
one coming from the volume and the other from the sur- 
face of the medium. These contributions are denoted by 
q^ and below. The total (that is, integral over the 
body volume) heat absorbed per unit time is given by 
the formula 



= <1 



( v \r)d 3 r + J> q( s \r)d 2 r , 



(5) 



where the first integral is evaluated over the body volume 
and the second over the surface. The quantity Q com- 
puted according to ((5|) is exactly the same as in the con- 
ventional theory. However, my calculations show that the 
volume contribution, oc Im(^e). I argue that in pas- 
sive media, the second law of thermodynamics requires 
that q( v ) > in contradiction with the inequality ((T|). 
In optically active media, it is possible to have < 
but the condition for negative refraction is then reversed 
and reads Im(^e) < 0. Thus I come to the conclusion 
that negative refraction is not possible in either passive 
or active media. 

The paper is organized as follows. In Section[2j I com- 
pute the volume contribution to the heating rate, q^ v \ 
for a monochromatic plane wave by two different meth- 
ods and obtain two different expressions. In Section 02 
I argue that the reason for this discrepancy is incorrect 
definition of the Poynting vector S. When the correct 
definition ([3]) is adopted, the two methods yield the same 
result. Also, in Section [3j is computed for general 
monochromatic fields (not necessarily plane waves). In 
Section [4] I compute the surface contribution to the heat- 
ing rate, q( s \ Also in this section, the zero-frequency 
limit is discussed. In Section [5] I give a detailed proof 
that the second law of thermodynamics requires that 
q(V) > o. In Section [6j I discuss compatibility of the ob- 
tained expressions for the heating rate with the causality 
principle. In Section [JJ I show that negative refraction 
is not possible even in anisotropic and nonlocal media. 
Finally, Sections [8] and [9] contain a discussion and a sum- 
mary of obtained results. 

2. Computation of the Heating Rate 

The heating rate q is defined as the energy absorbed and 
transformed into heat by a material per unit volume (or 
surface, if there is a surface contribution), per unit time. 
In the case of oscillating electromagnetic fields, this en- 
ergy must be averaged over time periods which are much 
larger than the characteristic period of oscillations. In 
this section, I use two different methods to compute q for 
a monochromatic plane wave propagating in a homoge- 
neous, isotropic medium characterized by scalar and lo- 
cal (but time-dispersive) functions e(ui) = e'(i-u) + ie"(uj) 
and = + ifi"(uj). 



A . First derivation of q 

First, we use the well-known formula for q which can be 
found in many standard textbooks, namely, 



(conv) 



1 

-in 



E 



<9D 

~dt 



H 



dB 

~dt 



(0) 



where (...) denotes time averaging and the quantities 
E, D and H, B are the electric field and displacement 
and the magnetic field and induction, respectively. The 
superscript "(conv)" has been used to indicate that ([6|) 
gives the conventional result for the heating rate. In the 
case of a monochromatic field of frequency u>, ([6|) can 
also be written as 
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(conv) _ I" // 



-in 



[e»<E 2 )+ M »(H 2 )] 



(T) 



Note that ©J® are quadratic in electromagnetic fields; 
correspondingly, E, D, H and B arc defined in these ex- 
pressions as real- valued quantities. 

Let us take one step further and evaluate (J7J) for a 
plane wave propagating in a homogeneous medium. We 
shall seek an expression for the heating rate which con- 
tains only the amplitude of the electric, but not of the 
magnetic, field. To this end, we write 



E = Re[Eoe^ k - : 



H = Rc[H e 



(k-] 



(8) 



where Eo and Ho arc complex field amplitudes. Time av- 
eraging yields (E 2 ) = (1/2)|E | 2 exp(— 2k" • r) and anal- 
ogously for the magnetic field. Here k" = Im(k) and k 
satisfies k • k = pie(ui/c) 2 . We now substitute the expres- 
sions for the time averages (E 2 ) and (H 2 ) in terms of 
the field amplitudes Eg and Hq into ((Jj to obtain 



(conv) 



8- 



[6»|Eo| 2 + m»|H | 2 ] e 



-2k"- 



(9) 



Further, we want to express |Ho| 2 in terms of |Eq | 2 - From 
the Maxwell equation cV x E = —d~B/dt and from Bo = 
/iHo. it follows that Ho = [c/w/x(w)]k x Eo. Therefore, 
|H | 2 = (c/o,| M |) 2 (k x Eo) • (k* x Eg) = (c/^H) 2 [(k • 
k*)(E • E^) - (k • E;;)(k* ■ E )]. The wave vector of a 
propagating (that is, not evanescent) wave can always be 
written as k = A:u, where u is a purely real unit vector 
such that u • u = 1 and k 2 = /je(uj/c) 2 is a complex 
scalar. In this case, k-k* = \[i(oj)e(uj)\(uj / c) 2 andk-Ep = 
k* - Eo = 0. The final expression for the heating rate then 
becomes 
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(conv) 



8^ 



ifMJ 



|Enl 2 e- 2k "' 



(10) 



Already at this point we can notice that the coefficient 
in the parentheses in the above formula appears to be 
somewhat strange. Indeed, if [i is purely imaginary (e.g., 
near a resonance), this coefficient becomes |e| + e". If, 
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in addition, |e'| 3> e", the heating rate becomes propor- 
tional to |e'|. 

We note that a propagating wave in an absorbing in- 
finite medium grows exponentially in the direction — k". 
To avoid the unbounded growth, one has to consider a 
half space z > into which an incident wave enters, 
e.g., from vacuum, and apply the condition z • k" > 0. 
However, a wave which is refracted from vacuum into 
an absorbing medium is necessarily evanescent. This fol- 
lows immediately from the fact that the projection of 
the wave vector on the plane z = must be contin- 
uous at the interface and, therefore, purely real (since 
it is real in vacuum). In the case of evanescent waves, 
the equalities k • k* = |/ie|(w/c) 2 and k • Eq = do not 
hold and the expression for g( conv ) becomes more compli- 
cated. This effect is not important for weakly absorbing 
media and it will not be discussed here. We only note 
that the expression (fT8|) which will be obtained below 
from the definition q = (J • E) applies to both running 
and evanescent waves and, in any case, differs from ([9]) 
or (Unj). 

B. Second derivation of q 

We now compute the same quantity as in the previous 
subsection but using a different, presumably equivalent, 
definition. Namely, we write 

?=(J-E), (11) 

where J is the total current in the medium induced by 
the propagating electromagnetic field. We again empha- 
size that this current is formed by the charged particles 
(bound and conduction electrons, ions, etc.) which make 
up the medium. Equation (|11| is simply the mathemat- 
ical formulation of the statement that, in a stationary 
state, the heating rate is equal to the (time-averaged) 
work that the electric field exerts on the medium per 
unit time per unit volume. 

We again consider a plane monochromatic wave with 
the electric field given by the first equation in ©. The 
current also has the form of a plane wave: 

J = Rc[J e' i(k - r - wt) ] . (12) 

Note that the above formula is valid only inside the 
medium volume. At the boundary, there is an additional 
surface current related to magnetization. This current 
and the corresponding contribution to the heating rate 
will be considered separately in Section [U We now focus 
on the volume contribution to the heating rate and de- 
note the corresponding quantity by q( v \ Time-averaging 
results in 

= Ir c (J • E*)e- 2k "' r . (13) 
To find Jo, we write the two curl Maxwell equations as 

cV xE= -dB/dt , cV x B = dV/dt + 4nJ , (14) 



where J is given by (J2J). Note that the above equations 
are equivalent to the usual macroscopic Maxwell equa- 
tions if we define the auxiliary fields D = E + 4irP and 
H = B 47rM and use ([2|) . By taking the time derivative 
of the second equation in (|14j) and substituting dB/dt 
from the first equation, we find that 

-4ndJ/dt = c 2 V x V x E + <9 2 E/di 2 . (15) 

At the next step, we substitute (|12[) and the first equa- 
tion in © into (|15p to obtain 

J = -(c 2 /47riw)[k x k x E + (w/c) 2 E ] . (16) 

We then use kxkx E = — (k-k)E and k k = iie{uj/c) 2 
(this holds for both propagating and evanescent waves) 
to further simplify the above expression for Jo, which 
becomes 

J = -^[/iMeH-l]E . (17) 

Am 

Upon substitution of the above expression into (fl~3|) , we 
arrive at 

(V) = ^W Im U^W^l e -2k"-r _ (18) 

I shall generalize this result to the case of monochromatic 
field E = Re[E w (r) exp(— iuit)] (not necessarily a plane 
wave) in Eq. (f3"0"|) below. 

The expression (fT5)l must be compared to (TTU1) . The 
respective formulae obviously differ. The reason for this 
discrepancy and the correct choice of the expression for 
q are discussed in Section [3] 

C. The two expressions for the heating rate and the 
constraints on e and fi that follow from them 

If we accept the conventional result for the heating rate 
as correct, the second law of thermodynamics requires 
that, in a passive medium, g( conv ) > 0, where g( conv ) is 
given for plane waves by Q or (fT0|) . If, however, we as- 
sume that the alternative expression (|18p is correct, then 
the second law requires that the volume contribution to 
the heating rate be positive (proof is given in Sec- 
tion [5]). From this, a different constraint on the possible 
values of e and /i is obtained. 

We note right away that for the conventional expres- 
sion (J9j) to be positive, it is not necessary that 

e" > , n" > , (19) 

although the above inequality is a sufficient condition. 
The sufficient and necessary condition is 

W + > . (20) 

Neither (fl"9")) nor PD|) prohibit negative refraction. We 
note that it can be argued (assuming © is correct) that 
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both inequalities e" > and f/' > must hold simulta- 
neously and independently to guarantee positivity of the 
heating rate [7,8]. Indeed, the condition ([20]) was derived 
for a plane wave and is, therefore, not the most general. 

The alternative expression (fTS)) imposes a different 
constraint on e and fi. For given by the expression 
(jTHj) to be positive, the sufficient and necessary condition 
is 

eV'+A">0. (21) 

Thus, e' and // can not be simultaneously negative while 
e" and fi" are positive. In particular, (|2"Tj) prohibits neg- 
ative refraction in the sense that if the wave number k 
satisfies k 2 = e/i(w/c) 2 , its real and imaginary parts have 
the same sign, independently of the choice of the square 
root branch. 

Finally, note the following interesting fact. If we put 
/j = I , formulae (jT7J|) and (fP8)) become identical. But if 
we put e = 1 (e.g., in a purely magnetic material), the 
two expressions still differ. This is suggestive of the fact 
that magnetic losses are not properly accounted for in 
one of these formulae. 

3. Correct Expressions for the Poynting Vector 
and the Heating Rate 

Both formulae (fTUj) and (fig)) do not obviously contradict 
any of the basic physical principles, such as the conser- 
vation laws. Therefore, we must choose the correct ex- 
pression for q on less fundamental grounds. To this end, 
we examine the origin of the two definitions ([6]) and (jlip . 

The expression (j6|) is obtained from q ~ — (V-S) where 
the Poynting vector S is given by Eq. Differentiation 
leads to 

q = — — (H • (V x E) - E • (V x H)) . (22) 

47T 

One then uses the macroscopic Maxwell equations to ex- 
press V x E and V x H in terms of the corresponding 
time derivatives to arrive at (|6|). In the stationary case, 
when there is no accumulation of electromagnetic energy 
anywhere inside the medium, the formula q = —(V-S) 
is undoubtedly correct. It is a mathematical expression 
of the statement that the total electromagnetic energy 
which enters into a small volume SV through its surface 
is entirely consumed to compensate for the irreversible 
(absorptive) losses in that volume. 

The alternative formula, (|11|) . is also a first-principles 
definition of the absorbed power per unit volume and 
appears to be unassailable. I emphasize again that the 
quantity J in is the total internal current produced 
by all charged particles that compose the material. This 
includes bound electrons, conductivity electrons (if such 
are present), ions in the case of plasmas, etc. 

So far, it appears that in either of the two approaches, 
the only formula that can be doubted is the definition of 
the Poynting vector ((4|) which was used to derive ((6|) . It 
should be noted that in standard textbook expositions, 



the form the Poynting vector is postulated rather than 
derived. Thus, for example, Schwinger et al. (in Classical 
Electrodynamics [9, $7.1]) consider the identity 

^V. (E ,H, + ±( E .§ + H.f)=0 ,23, 

which is trivially obtainable from the macroscopic 
Maxwell equations in the absence of external currents. 
Then Schwinger et al. write: "Our aim is to write this 
result as a local energy conservation law. We immedi- 
ately identify, from the divergence term, the energy flux 
or Poynting vector S to be S = (c/47r)E x H." The ar- 
gument is, however, mathematically flawed. Indeed, one 
can take any scalar function /(r) ^ whose integral over 
the body volume is zero and write it as a divergence of a 
vector field, /(r) = V-F(r), where F(r) vanishes outside 
of the body. One then can add (c/47r)V • F to the first 
term in the left-hand side of ([23]) and subtract it from 
the second term, and the identity will still hold. Accord- 
ing to the logic of Rcfs. 9, one then has to define the 
Poynting vector as S = (c/4tt)[E x H + F]. Note that 
the field F does not need to be solenoidal, so that not 
only the definition of S is changed, but also of its diver- 
gence. This ambiguity in the conventional definition of S 
leads to a substantial strain. To quote Schwinger again, 
"... More intractable is the identification of the last term 
in (|2"5|) ." I argue that such identification is, indeed, in- 
tractable because the term in question has no physical 
meaning. 

A somewhat different approach to deriving the con- 
ventional expression for S is adopted by Landau and 
Lifshitz in Electrodynamics of Continuous Medium [10, 
$80]. First, it is shown that Eq. (|4]) is valid in non- 
magnetic media where H = B. Then Landau and Lif- 
shitz argue that the normal component of S should be 
continuous when a wave crosses an interface between two 
media. Since the tangential components of both E and 
H arc continuous, the normal component of S defined 
by ([4]) is continuous as well. Therefore, (J4| should be 
valid in any media, including those with dispersion and 
a magnetic response. 

I do not dispute here that the tangential components 
of E and H are continuous as long as there is no surface 
current at the interface which is formed by charges which 
are external to the medium. Note that sometimes such 
currents are referred to as "free currents" (even though 
they do not include the current of free electrons in the 
case of conductors). However, I claim that continuity of 
the normal component of S is an incorrect boundary con- 
dition for interfaces that separate two media at least one 
of which is magnetic. Indeed, it is known [10, $29] that 
nonzero magnetization creates surface currents which arc 
restricted to a very thin layer near the medium boundary. 
These currents are formed by the charges of the medium 
and, therefore, do not cause the tangential component of 
H to become discontinuous as follows immediately from 
the equation cVxH = dD/dt. When a wave crosses an 
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interface in which such surface current is flowing, a finite 
fraction of its energy is lost to the (positive or negative) 
work exerted by the electric field on the surface current. 
In this case, the normal component of S experiences a 
discontinuity. The role of the surface currents and their 
input to the heating rate is discussed in Section [4] below. 

To obtain the correct expression for the Poynting vec- 
tor, we start with the microscopic electric and mag- 
netic fields, e and h. The spatial averages of these fields 
are [10, $1 and $29] e = E and h = B. Here the bar 
denotes spatial averaging over physically small volumes. 
Further, we write e = E + <5e and h = B + Sh, where 
5e and Sh are the fluctuating parts of the fields. The 
microscopic expression for the Poynting vector is 



7 (V) 



-v-(s> 



-Re[-V-(E;xB 

c 



ReK.(VxB„)-B„.(VxE:)] 



-Re 
-Re 



8- 

c 

c 

8tt 
^|E^| 2 
8tt 



E* • (V x M^)HJ 



-B„, • B ,., 



-IU) 



Im[/z(w)e(w)] 



(29) 



We thus have derived the following formula for q^: 



s = — e x h 

47T 



We now average the above expression as follows: 



(24) 



S =s = 



-in 



[ExB + <5exB + Ex<5h + <Sex<5h] 



(25) 



By definition, Se = Sh = 0. The term Se x Sh is 
quadratic in field fluctuations and can be omitted as 
small. It should be also noted that in materials which 
are random but isotropic on average, this term is identi- 
cally zero by symmetry. We thus arrive at the expression 
(J3j> for the Poynting vector. We note that in vacuum, the 
Poynting vector is expressed in terms of the electric and 
magnetic fields. There is no conceivable physical reason 
why this should change if the field propagates through a 
material medium. But the average value of the magnetic 
field in the medium is B, not H. Despite the fact that 
H is commonly called the " magnetic field" , it is actually 
an auxiliary quantity. 

Let us adopt the definition ([3]) for the Poynting vector 
and compute from q' y ) = — (V ■ S) for a general 
monochromatic fields of the form 



E = Re[E w (r)e^ t ] , D = Re [D w (r)e^ 4 ] , (26) 
H = Re [H w (r)e- Mt ] , B = Re [B w (r)e~ iwt ] , (27) 

where D w = e(w)E w , B w = ^(w)H w , cV x E w = iuB^ 
and cV x H w = — iusD^. At the moment, we do not 
consider the heating rate at the surface where S has a 
discontinuity. We then have: 



(s) = x B a 



(28) 



q[v) = -|^ Im[MHeM] 



(30) 



If we set E w = E exp(zk • r), the above expression co- 
incides with formula (|18p which was derived previously 
from the definition = (J • E) for the case of a plane 
wave of the form ©. 

Next, consider the definition qW) = (J • E). We have 
already used this definition to compute for a plane 
wave with the result given by Eq. (fl~8|) . Now we repeat 
the calculation for more general monochromatic fields 
((26|) . ((27l) . The current (except at the medium surface) 
can also be written in a similar form, namely, 



J = Re [J w (r)e- Mt ] 



where 



We then write 



4tt 



(cV xB.I iuj^ u 



(31) 



(32) 



,(V) 



(J • E) = -Rc(J w • E*) 

-^Re [iuj-E^ ■E:+c(VxB u ). e; 

07T 



-Re[^(w)(V x n u 



E* 



87T 
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8^r 
wlEJ 2 



-Re 



- M (w)e(w)E w • E; 



8tt 



Im[/j(w)e(w)] 



(33) 



The result of this calculation coincides with (|30l) . 

Thus, we can conclude that if S is defined by ([3]), the 
two definition of the heating rate, q^ = — (V • S) and 
q( v ) = (J-E) are equivalent and we have the statement of 
local energy conservation which, in the stationary case, 
reads 



and 



(J • E) + (V ■ S) = 



(34) 
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4. Heating Rate at the Surface and the Total 
Absorbed Heat 

The derivation of the heating rate form the formula (fTTj) 
was so far restricted to points inside the medium. In this 
section, the additional surface term q( s 1 is derived. 

First, consider the definition (fTTj) . In the monochro- 
matic case, the current in this formula is given by (|32[) 
where B w = /i(w)H a) . In the sequence of equalities (|33[) . 
I have, at one point, replaced the term V x B w by 
/x(o;)V x H w . This operation is only valid inside the 
medium volume. Close to the surface, we must write 

VxB u = Vx mO)H w = /u(r)V x H w + [V/i(r)] x H u , 

, (35) 

where the dependence of fj, on position has been indicated 
explicitly In the case of macroscopically homogeneous 
media, V/i(r) = everywhere except at the surface, 
where ^i(r) experiences a discontinuity. If we restrict at- 
tention to points r which arc on the surface, evaluation of 
V^t(r) results in the additional surface current [10, $29] 

J(f) = cn x M w , (36) 

where n is the outward unit normal to the boundary 
at the point r. Note that the definition of the surface 
current (|36|) does not contain a spatial delta-function and 
that has different physical units then the volume 
current J. 

We now find the surface contribution to the heating 
rate as 

= 1 Re ) . E *) . (37) 
A straightforward derivation results in 

= ^ Rc K 1 - M)( n x H -) ' E -l • (38) 

It is also possible to start from the definition q = — V • S, 
take into account the fact that the normal component 
of S defined by ^ experiences a discontinuity at the 
medium boundary, and arrive, in a straightforward man- 
ner at 

9 (S) = ^Rc[(l- M )(H w xE:).h] . (39) 

Since a • (b x c) = b • (c x a), the two expressions ([38)1 
and (f39j) are identical. 

The total heat absorbed by the body, Q, is given by 
Eq. (0) • It is easy to see that this quantity is the same as 
in the conventional theory. Indeed, Q can be computed 
by integrating the energy flux through any surface en- 
closing the body. Such surface can be drawn in free space 
where the conventional expression for the Poynting vec- 
tor (0| and the expression derived in this paper © co- 
incide. Therefore, the proposed change in the form of the 
Poynting vector and of the heating rate does not affect 
any of the previously established results for differential 
or integral cross sections, such as the Mie formulae for 



extinction, absorption and scattering cross sections of 
spheres. 

The surface contribution to the heating rate derived 
in this section requires several additional comments. 
The obvious distinction between the surface term q^ 
([39]) and the volume term q^ (f30|) is that the volume 
term is proportional to the frequency to while the sur- 
face term is not. Of course, it is incorrect to say that 
q( v } always vanishes in the zero-frequency limit because 
lim UJ ^o[ w e(w)/i(a;)] = Airia^uj = 0), where a is the 
static conductivity of the material [11]. However, it ap- 
pears that the surface term does not vanish in the zero- 
frequency limit even if we formally set a = 0. This pos- 
sibility is worrisome and is discussed below. 

First, this paper is primarily concerned with the high- 
frequency superficial magnetism which originates due to 
the loop-like conductivity currents flowing in elementary 
cells of composite materials. The magnetic susceptibil- 
ity of such composites identically vanishes in the zero- 
frequency limit [12,13]. Moreover, the formula (fTTj) in 
which the current is given by © assumes that all cur- 
rents obey the classical laws of motion. This may not be 
the case when magnetization is caused by spin aligning, 
as in the cases of ferro- and para-magnetism. However, 
even in the case of ferromagnetism, the current cV x M is 
a macroscopic quantity. According to the Ehrcnfcst the- 
orem, all macroscopic quantities obey the classical laws 
of motion. Nevertheless, it should be acknowledged that 
some of the phenomena associated with ferromagnetism, 
such as the hysteresis, are clearly outside of the theoret- 
ical frame of the classical electrodynamics of continuous 
media. Such effects can be accounted for phcnomcnolog- 
ically but not in a fully self-consistent way. 

Thus, the zero-frequency limit might not be the proper 
test for the theory developed in this paper. I will, how- 
ever, argue that it is possible to apply this theory to the 
zero- frequency limit without obtaining unphysical effects 
or contradictions. To this end, I consider below two sim- 
ple examples. 

The first example is a straight ferromagnetic or para- 
magnetic cylindrical wire of radius a, conductivity a and 
permeability fi carrying a current of uniform density J 
directed along the axis of the wire. I disregard here the 
Hall effect that results in a non-uniform current distribu- 
tion over the wire cross section [14]. The quantities fj, and 
a are purely real at zero frequency. Then my theory pre- 
dicts that the volume will be heated at the rate (per 
length L of the wire) 

Q^/L = ira 2 anE 2 (40) 
and the surface will be heated or cooled at the rate 

Q (s) /L = ira 2 a{l - n)E 2 . (41) 

In a ferro- and paramagnetic materials, fi > 1 (for the 
case of ferromagnetic, the nonlinearity of the magneti- 
zation curve and the magnetic memory of the material 
must be taken into account, which is not a trivial task), 



G 



so that the surface term is negative. But the total heat 
produced in the system is the sum of both contributions, 
namely, 

Q/L = Q (y > /L + Q {S) /L = ira 2 aE 2 , (42) 

which is the Joule's law. In a steady state, the overall flux 
of thermal energy through the wire surface, which is the 
experimentally measurable quantity (e.g., in a calorime- 
ter) is given by Q/L, in agreement with the Joule's law. 

The second example is a magnetic object placed in 
crossed external electric and magnetic fields. For sim- 
plicity, consider a long cylinder uniformly magnetized 
along its axis. The magnetization will create loop-like 
surface currents that flow around the cylinder axis. We 
now place the cylinder in an external electric field which 
is perpendicular to the cylinder axis (this example has 
been previously considered by Pcrshan [15]). If the cylin- 
der is conducting (as are most fcrromagnets) , the tan- 
gential component of the electric field at the cylinder 
surface, as well as the electric field inside the cylinder, 
vanish and we obtain = q( v *> = 0, as expected. A 
somewhat more complicated situation arises if we for- 
mally set (j = and [i — 1 ^ 0. The volume term q( v > is 
still zero in this case, but the surface term q^ s \<p) may 
become locally nonzero (here p, z,tp are the cylindrical 
coordinates). Even though it can be easily seen that 

f q {S \^ = 0, (43) 
Jo 

we still expect no local heating or cooling of the surface 
in the static equilibrium. The contradiction is resolved 
by noting that the state of the cylinder described above 
can not be its true state of equilibrium and that the 
initial assumption a = and ^ was unphysi- 

cal. This assumption contradicts the mechanical equilib- 
rium of charges that make up the circular magnetization 
currents. Classically, these charges rotate with constant 
angular velocity around the cylinder axis due to a phc- 
nomenological radial force. It is, however, not possible to 
introduce a phenomenological restoring tangential force 
such as the harmonic restoring force in the Lorentz model 
of dielectrics. Indeed, such restoring tangential force will 
preclude the magnetization current from flowing in the 
first place. Consequently, imposition of an external tan- 
gential force (due to the external electric field) in the 
absence of a tangential restoring force is bound to break 
the equilibrium of the system. Specifically, the external 
electric field will cause electric charge to accumulate on 
the cylinder surface until the tangential component of 
the electric field is completely nullified. The resultant 
state will be the true static equilibrium of the system. 
The conclusion is that magnetized objects can not have 
identically zero conductivity. Of course, the value of a 
can be small, but so is usually the value of fi — 1. An- 
other important consideration is that, in addition to con- 
ductivity, magnetics also have some dielectric response 
whose effect is to diminish the tangential electric field at 



the body surface. 

The example considered above suggests that the con- 
ventional definition of the Poynting vector (j4| is erro- 
neous because it predicts existence of an equilibrium 
state which contradicts mechanical stability of the sys- 
tem. 

5. Thermodynamic Considerations and Impos- 
sibility of Negative Refraction 

Many authors believe that the unique properties of the 
negative refraction materials originate from the fact that 
the phase velocity and the Poynting vector in such media 
are oppositely directed. This property is sometimes re- 
ferred to as "backward propagation" . For example, to 
quote Marques at al. [16, $1.2], "...most of the sur- 
prising unique electromagnetic properties of these me- 
dia arise from this backward propagation property." If 
the expression ([3]) for the Poynting vector is correct, as 
I argue in this paper, then the phase velocity and the 
Poynting vector always point in the same direction and 
"backward propagation" is impossible, at least in elec- 
tromagnetically homogeneous media. I will, however, ap- 
ply more fundamental thermodynamic considerations to 
show that the inequality §Q is physically prohibited, re- 
gardless of whether it results in those "surprising unique 
effects" or not. 

To this end, it is instructive to introduce the " accessi- 
ble heat". This is the heat (either positive or negative) 
which can be transferred from the body to a heat reser- 
voir on a time scale which is short compared to the time 
scales associated with heat diffusion in the body. Obvi- 
ously, this is the heat generated at the surface. Let 

Q(S) = j q (S) {r)d 2 r = Q {S) _ Q (S) (44) 

Here Q + is obtained by integration over the surface 
areas where q^ s \r) is positive and is obtained by 
integration over the surface areas where q^ s ' (r) is neg- 
ative. Let us further assume that the material exhibits 
negative refraction and q^^r) is negative, so the heat 
generated in the volume, 

q(v) = f q (v) {r)d 3 r (45) 
Jv 

is also negative. Then we have 

Q = -\Q (V) \+Q { + 1 -Q { - ] (46) 

or, cquivalently, 

Q { + ] =Q + \Q (V) \ + Q (S) >Q ■ (47) 

Thus, the positive accessible heat is greater than the 
total heat absorbed in the body. I will now demonstrate 
that this contradicts the Carnot theorem and, moreover, 
can be used to create a perpetuum mobile of the second 
kind. To see that this is, indeed the case, consider the 
cyclic process shown in the Fig. Q] 
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Fig. 1. A cyclic process involving a negative refractive 
index material that violates the Carnot theorem. The 
black oval represents a negative refraction medium and 
the white oval is an ideal Carnot engine. 



In this cycle, the following events happen: (1) A 
negative-refraction sample represented by the black oval 
(referred to as the "body" below) at the initial tempera- 
ture Th is irradiated for a period of time At which is 
short compared to the time scales associated with heat 
diffusion in the body, yet long compared to the electro- 
magnetic oscillations period, so that radiation is almost 
monochromatic. The body absorbs the energy 8 = QAt 
from the radiation field. (2) The body is brought in con- 
tact with a heat reservoir at the temperature Th which 
has very high heat conductivity; the amount of ther- 
mal energy 8^ = At generated at the body's sur- 
face is transferred adiabatically to this reservoir. (3) The 
body is disconnected from the reservoir and heat diffu- 
sion takes place in the body until the equilibrium tem- 
perature T 1 < Th is reached. (4) An ideal Carnot engine 
is operated for one cycle between the hot reservoir and 

a colder reservoir whose temperature is Tl < Th- The 

( s) 

Carnot engine absorbs the heat 8\ from the hot reser- 

Vc9+^ and rejects some 



voir, makes useful work A 



amount of heat 8' = 8 



(S) 



A to the cold reservoir. Here 



Vc = 1 - T L /T H 



(48) 



is the efficiency of an ideal Carnot engine. (5) The body is 
again brought in contact with the hot reservoir; now the 
heat | flows back from the hot reservoir to the body. 
In the end of this process, the body has the temperature 
Th- (6) We disconnect the body from the hot reservoir. 
Now the cycle is complete and the system has returned to 
its original state. Note that steps (I) and (2), as well as 
(3) and (4) can be combined, i.e., occur simultaneously. 
The net effect of the above thermodynamic transfor- 



mation is the following: The electromagnetic field has 
done the work 8 on the body which was immediately 
dissipated into heat 8] we then converted this heat into 
the useful work A. The overall efficiency of this process 



A 

V = T = VC- 



3 (S) 



\e<y)\ 



i - \o<y)\/o. 



(S) 



(49) 

in violation of the Carnot theorem. Moreover, we can 
operate a perpetuum mobile of the second kind if A > 8 
or, equivalently, r] > 1. This is achieved if Tl/Th < 
|0(v)|/0(S). xhere 

is no physical reason why this condi- 
tion can not be met. In particular, it can be met quite 
easily in the case of low-loss negative refraction materi- 
als such that \8 ( - v '>\/8^ ) = 1-6 where 8 < 1. Then even 
relatively small temperature difference Th — Tl would 
be sufficient to extract more energy from the heat reser- 
voir than was absorbed from the electromagnetic field. 
Note that in order to obtain the contradiction, it is es- 

(S) 

sential that the "accessible" heat 8\_ be larger than the 
total absorbed heat 8. This is always the case for nega- 
tive refraction materials, as was shown above. Therefore, 
I conclude that negative refraction is impossible. 

It is also possible to use negative refraction to con- 
struct a refrigeration cycle in violation of the Carnot 
theorem. 

6. Heating Rate and Causality 

As is well known, the causality principle requires that 
e(u>) and n(u>) have no singularities in the upper half- 
plane when viewed as functions of the complex vari- 
able u>. In the lower complex half-plane, these functions 
may have singularities. We will assume now that all 
such singularities are simple poles. If, in addition, we 
account for the symmetry property e(— u>) = e*(w) and 
fj,(— to) = /j,*(u>), we can write e and \i in the most general 
form as 



e(o)) = I + 47rx e H , fj,(w) = 1 + 4nx m (u) , (50) 
where 



k 



AM 



a\ — (jJ 2 — iakUJ 



b\ — uj 2 — ifikUJ 



(51) 
(52) 



Here \ e and Xm ar e the electric and magnetic suscepti- 
bilities, respectively, afc, bk, ctk and /3k are coefficients and 
/fe(w) and (?fe(w) are analytical functions of the frequency 
which have no singularities and satisfy fk(—u>) = //*(w) 
and analogously for gk- The representation ((50|) -([52 |) is 
customary in the theory of dispersion. In the electric 
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case, the functions fk are usually positive constants in- 
terpreted as oscillator strengths. In the magnetic case, 
the typical form of <?fc(w) is gk{^>) = c^w 2 where the 
coefficients Ck can be negative, i.e., in the case of dia- 
magnetics. Terms proportional to ui 2n with 77 > I in the 
Taylor expansion of /fc(o>) and gk{u) arc not usually con- 
sidered because of the physical requirement that both e 
and (j, are bounded when uj — > 00. The physical interpre- 
tation of the remaining constants appearing in (|5ip , (|52| 
is as follows: ak and bk are the resonance frequencies and 
ctk and Pk are the respective relaxation constants. 

We wish to examine whether the representation 
(f50|l - (|52|) is compatible with the inequality q(u>) oc 
lm[e(ui) h(uj)] > which must hold for all positive fre- 
quencies. Obviously, the latter imposes some constraints 
on the coefficients appearing in formulae (|5"Tj) - (f5"2")) . How- 
ever, it is easy to show that coefficients that satisfy 
q(ui) > for all u > do exist. For instance, a sufficient 
condition for q > is obtained when one of the suscep- 
tibilities is significantly smaller than the other. Let us 
write 



Im(e/x) = 47rlm(xe + Xm) + (47r) 2 Im(x e Xm) • (53) 



Since (4n) 2 lm(xeXm) > ~{^) 2 \XeXm\ and lm(x e + 
Xm) > Ini(x e ), a sufficient condition for Im(/ie) > is 
\Xm\ < Xe/^lXel- This inequality can always be satis- 
fied in the whole frequency range for sufficiently small co- 
efficient Ck (assuming <7fc(u;) = c^u 2 ). Analogously, q > 

if IXe| < X'L/^lXml- 

Thus, we have obtained two sufficient conditions for 
q > 0. However, these conditions are not necessary. To 
derive a condition which is both sufficient and neces- 
sary, consider the special case of a single electric and 
single magnetic resonance with positive and frequency- 
independent oscillator strengths, /f and g 2 . That is, as- 
sume that e and [i are given by 



e(w) = 1 + 4tt 



fl(w) = 1 + 47T 



UJ 2 — UJ 2 — IJeUJ 



<x 2 rn - J 2 - ij m u> 



(54) 
(55) 



where all coefficients arc positive. Then a straightforward 
calculation shows that 



A e (m)A m (uj) 

A-KLJ 



Im(fie) — auj 4 + but 2 + c , (56) 



where 



a 
b 



fe Te "I" fmlra ? 

/£ 7m (7e ~ 2^ e 2 ) + fl^l - 2u 2 m ) 

-47T/ e 2 /^(7e+7m), 
47r /e /m(7mW e 2 + Je^m) 



(57) 



(58) 



(59) 



and A e (u)),A m (w) are positive factors defined by 

A e (L0) = (lo 2 - lo 2 ) 2 + ( le uj) 2 , (60) 
A m (uj) = {oj 2 m - to 2 ) 2 + ( 7m w) 2 . (61) 



The right-hand side of (|56j) is a quadratic polynomial 
in lo 2 with a positive free term c. If we assume that 
relaxation constants are small so that j 2 < 2lo 2 and 
j m < 2w 2 j, the coefficient b is negative. In this case, the 
necessary and sufficient condition that the polynomial is 
positive is V = b 2 — Aac < 0. A tedious but straightfor- 
ward calculation yields the following expression for the 
discriminant T>: 



V 

7e7m 



lelm [fmile ~ ^e) + f e { lm ~ 4u, m )] 



2 2 
'7e w m 



+Vtf m [llll ~ 2(7™^ 2 

-2(^ + ^ + 4^1] 
-8vr/ e 2 /^ {/* [ 7e 2 + 7e7m + 2(w 2 n - 

+ fe [7,i + 7rn7e + 2K 2 
+ (^) 2 fef m (le+lm) 2 ■ 



<)] 



,2 



^m)] } 

(62) 



Assuming that b < 0, the sufficient and necessary con- 
dition for Im(/ie) > is that the above expression is 
negative. This, of course, leads to a very complicated in- 
equality However, if the relaxation constants are small 
compared to all other physical scales of the problem, this 
inequality can be simplified and reads 



(c 2 -^) 2 +4 7 r(^-c 2 )(/ 2 i -/ 2 ) 



> 4tt 2 / 2 / 1 



P 2 ( le + 7 m ) 2 
7e7m 



(63) 



To illustrate how inequality (|63[) works, I have plot- 
ted Im(/ie) as a function of u> for two different sets of 
parameters (see Fig. [5] caption for details). As can be 
seen, if the parameters satisfy (|63|) . Im(/ie) > in the 
whole frequency range (Fig. [5^). If, however, we use a set 
of parameters that does not satisfy (|63[) . there appears a 
frequency range in which Im(/je) < (Fig. [2b)- As I have 
argued above, q oc lm(pe). Therefore, the set of param- 
eters used to compute the curve shown in Fig. [2Jd leads 
to negative heating rate in the frequency range indicated 
by the horizontal arrow. These values of parameters vi- 
olate the second law of thermodynamics and, therefore, 
can not be realized in any material, either natural or 
artificial. 

Thus, we have seen that the expression for the heating 
rate derived in this paper does not contradict causality. 
However, it imposes constraints on the possible values of 
constants in the dispersion formulae (fSTj) . (f52"|) . Generally, 
these constraints are very complicated mathematically. 
In the simplest case of one magnetic and one electric 
resonance (e and \x given by formulae (|54[) and (1551) ). 
the condition is that the discriminant (|62[) is negative. 
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In the limit of small relaxation constants 7 e and j m , this 
condition can be approximated by the much more simple 
inequality (l63l) . 

We finally note that the condition |63|) is only suffi- 
cient but not necessary if the relaxation is so strong that 
the coefficient b given by formula (|59|) is positive. Also, 
in the case when one of the resonance frequencies is zero, 
as is the case for electric permittivity of conductors, in- 
equality (|63j) may become an inaccurate approximation 
of the more general inequality T> < where T> is given 
by (52). 
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Fig. 2. Plots of Im[/j(w)£(w)] where e and fi are given 
by (|54|) . (|55|) for different sets of parameters. Plot (a): 

LO m = W , LO e = 4w , fm = 0.5w , fe = ^0, 7m = 0.01w 

and 7 e = O.lwo; wo is an arbitrary frequency scale. The 
parameters satisfy inequality f|63[) . Plot (b): same pa- 
rameters as in plot (a) but uj e = 2ujq. With this change, 
(|63p is no longer satisfied in the frequency range denoted 
by the horizontal arrow. Negative values of Im[/i(w)e(w)] 
are not shown in the plots due to the use of logarithmic 
scale. 



7. Heating Rate in Anisotropic and Nonlocal 
Media 

The general case of a medium with magnetic and elec- 
tric anisotropy and nonlocality is quite complicated. The 
wave vector k of a plane wave propagating in such a 
medium can be found from the following condition: 



det 



e 'kx/i 1 k x 



c 







(64) 



where e = e(w, k) and p, — k) are k-dependent ten- 
sors. 



The dispersion relation (|64[1 is simplified for the case 
of propagating waves. A propagating (as opposed to 
an evanescent) wave is characterized by a wave vector 
k = fcu where k is a complex scalar and u is a purely 
real unit vector such that u • u = 1. Thus, a propagating 
wave can, in principle, experience spatial decay. The im- 
portant point is that, in the propagating case, the wave 
vector is completely characterized by a direction in space 
(the unit vector u) and by a single scalar k. We can uti- 
lize this property to rewrite (|64p as 



det 



-k z T + ( - 



= 0, 



(65) 



where f = — e _1 u x /"r'ux. In general, the 3x3 tensor 
T is symmetric but not Hcrmitian. Therefore, its eigen- 
vectors and eigenvalues, denoted here by Vj and 1/tj, 
are complex. For each direction u, the wave number of 
a propagating wave is determined from one of the equa- 
tions k 2 = Tj(u>/c) 2 while the polarization of the j-th 
mode is given by E = avj, a being an arbitrary com- 
plex constant. 

Interestingly, it is possible to make a statement about 
the restrictions that are imposed by the condition q( v ^ > 
on the wave number k without explicitly solving the 
dispersion equations. We note that the formula q = 
(J • E) is valid in any elcctromagnctically homogeneous 
media. We then consider monochromatic, propagating 
plane wave with the wave vector k, so that the fields are 
of the form 



E = Re 
B = Re 
J = Re 



Eoe i(k-r-o<t) 

Boe i(k-r- wt ) 
j i(k.r-wt) 



and 



47rJo = iwEo + ick x Bo 
We then obtain 



7 (V) 



-2k"r 



8- 



-Im[(kxB )-E*] 



(66) 
(67) 
(68) 

(69) 
(70) 



We now use kxBo= (c/w)k xkxEo and the identity 
a x b x c = b(a • c) c(a • b) to arrive at the following 
result: 



,-2k". 



$7r(w/c) : 



-Im [|Eo| 2 (k.k)-(k-Eo)(k-E5)] 



(71) 

Note that we have not used any constitutive relations in 
the derivation of (|7Tj) . Also note that the wave vector k 
must satisfy the dispersion relation . If k = ku, it is 
always possible to write (k • E )(k • Eg) = cos 2 6>|E | 2 , 
where 9 is a purely real angle. From this, we obtain the 
final expression for q: 
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(v) _ ^|Eo| 2 e - 2k "- r S in 2 gIm(fc 2 ) 
q ~ Stt {uj/cY ■ {U) 

The phenomenon of negative refraction requires that 
the direction in which a wave exponentially decays due 
to absorption in the medium is opposite to its phase 
velocity. Mathematically, this means that the real and 
imaginary parts of the complex wave number must have 
opposite signs. But for this to be true, it is required that 
Im(fc 2 ) < 0. However, Eq. (J72J) implies that Im(fc 2 ) > 
0. Thus, negative refraction is not physically attainable 
even in anisotropic and nonlocal media. 

Finally, we note that one can formally choose the 
polarization and the wave number in such a way that 
sin 2 9 = so that the medium does not absorb radiation. 
However, it is easy to see that waves with sin 2 9 = do 
not satisfy the dispersion relation (|64|) . 

8. Discussion 

In this section, I address certain anticipated objections 
to the theory developed in this paper, as well as discuss 
some of its limitations. 

The first and the most obvious objection is that there 
have been a number of works which claim experimental 
demonstration of negative refraction in clcctromagncti- 
cally homogeneous materials. Such experiments can be 
classified into two kinds. The direct-kind experiments 
measure the deflection of a beam passing through an ex- 
perimental sample made of a subwavclcngth-structured 
"metamaterial" . However, in most experiments of this 
kind, the linear size of the smallest metamaterial ele- 
ment, £ is not much smaller than even the vacuum wave- 
length A. Strictly speaking, I should be compared to 
the wavelength inside the material. Further, the more 
physically relevant parameter is k£ = 2ir£/\. In typi- 
cal experimental demonstrations of negative refraction, 
this parameter is of the order of unity. Under these cir- 
cumstances, interpretation of the experimental results 
in terms of the bulk constants e and fi is problematic. 
There are also experiments in which the negative re- 
fraction is measured indirectly by means of measuring 
the transmission and reflection coefficients t and r of a 
subwavelength-structured thin film. Here, as in the case 
of direct-kind experiments, it is very difficult to achieve 
k£ <C 1. Additionally, the indirect-kind experiments rely 
on phase measurements of high-frequency electromag- 
netic fields and on a analytical procedure of extract- 
ing e and fi from the measurements of t and r. Both 
of these tasks are notoriously difficult and have recently 
been subject to some controversy [17-19]. 

The second objection is based on the factually incor- 
rect, yet widespread belief that the magnetic field can, 
under certain circumstances, do work, i.e., on magnetic 
moments. Theory developed in this paper is based on the 
premise that only electric field can do work. The question 
of whether the magnetic force can do work is simultane- 
ously simple and complicated. Of course, it immediately 



follows from the expression for the Lorentz force that the 
magnetic force does no work on a moving charged par- 
ticle. Yet, apart from this simple observation which can 
be found in most textbooks on classical electrodynamics, 
there has been almost no serious discussion of this ques- 
tion in scientific literature. At the same time, situations 
in which the magnetic force is apparently doing work arc 
quite abundant. Recently, the question was addressed 
in a mathematically rigorous way by Dcisslcr [20]. In 
this reference, it is shown that the magnetic force does 
no work on a classical magnetic moment under any cir- 
cumstances. It is further shown that the magnetic force 
does no work on an atom where magnetization is due to 
orbital angular momentum. Finally, Deissler shows that 
there is a fully self-consistent description of the quantum 
spin in which the magnetic field does no work either. 

Third objection is based on the belief that compos- 
ite media can be assigned effective medium parameters 
which describe (approximately) some phenomena associ- 
ated with wave propagation through such media but not 
the others. I believe that such contention was expressed, 
for example, by Simovski [19], although implicitly. My re- 
ply to this is that most experimentally measurable quan- 
tities, such as the intensity, are bilinear in the electric 
and magnetic fields. Therefore, any useful homogeniza- 
tion model must correctly predict such quadratic combi- 
nations, including the Poynting vector and the heating 
rate. 

Fourth objection is that the "magnetic" current cV x 
M is somehow different in its physical properties from 
the "electric" current dP/dt and, therefore, obeys dif- 
ferent laws of motion. Of course, in the case of meta- 
materials, both currents have exactly the same physi- 
cal origin. But even in the most general case, both cur- 
rents are macroscopic and there is no valid physical ba- 
sis to apply different laws of motion to them. It is also 
not possible to do so mathematically Assume that we 
know a vector field J(r). Assume also that we know 
that J(r) = J e (r) + J m (r) where J e (r) = dP(r)/dt and 
J m (r) = cV x M(r). Is it possible to find uniquely J e 
and J r „ if we know J (but not P or M)? It is known that 
J TO is solenoidal. If it were also known that J e is irrota- 
tional, we would be able to use the Hclmholtz theorem 
to uniquely decompose J into the irrotational and the 
solenoidal parts corresponding to J e and J m . But the 
only instance when J e is irrotational is the static case. 
Therefore, beyond strict statics, there is no unique way 
to disentangle the term cV x M from the total current 
J. 

The fifth set of objections is related to the zero- 
frequency limit. This limit is discussed in detail in Sec- 
tion [U Here I would like to reiterate the following. The 
theory developed in this paper is based on the fundamen- 
tal assumption that all currents obey the same classical 
laws of motion. Although the author sees no physical rea- 
son for this assumption to be untrue even in the zero fre- 
quency limit, reasonable caution must be exercised when 
applying the results to magnetization caused by quan- 
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turn spins. It is theoretically possible that the magnetic 
susceptibility which is due to spin alignment has different 
physical and mathematical properties when compared 
to the magnetic susceptibility which is due to classical 
eddy currents. In this case, the total permeability must 

be Written as \l - 1 = (^classical - 1) + (^quantum - 1), 

where //quantum is the contribution to the total perme- 
ability due to quantum effects which arc manifest only 
at low frequencies, and /i should be replaced by /^classical 
in the expressions for the heating rate derived in this 
paper. 

Finally, the essential requirement for applicability of 
the results derived in this paper is that the medium 
is electromagnetically homogeneous or can be effec- 
tively homogenized. Mathematically, this means that the 
medium must support running plane waves as its electro- 
magnetic modes. This condition is not satisfied in pho- 
tonic crystals and similar structures which support prop- 
agating modes in the form of Block waves. 

9. Conclusions 

The article has the following conclusions: (i) The correct 
definition of the Poynting vector in magnetic media is 
(J3j) - If this definition is used, then the local energy con- 
servation law in the form (|34|) holds, where J is given by 
fl2J. (ii) The heating rate q^ v \ defined as the amount of 
energy converted to heat per unit time per unit volume, 
is proportional to the factor Im(fie) inside the volume 
occupied by the material] heating rate at the surface is 
given by Eqs. ([38]) or ([39]) in Section H (iii) It follows 
from the detailed thermodynamic considerations of Sec- 
tion [5] that negative refraction contradicts the second 
law of thermodynamics. This statement holds for active 
or passive media and in the presence of anisotropy and 
spatial dispersion, as is shown in Section [7l 
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